0 Ju n 20 00 Fourth Order Algorithms for Solving the Multivariable
نویسنده
چکیده
We develop a fourth order simulation algorithm for solving the stochastic Langevin equation. The method consists of identifying solvable operators in the Fokker-Planck equation, factorizing the evolution operator for small time steps to fourth order and implementing the factorization process numerically. A key contribution of this work is to show how certain double commutators in the factorization process can be simulated in practice. The method is general, applicable to the multivariable case, and systematic, with known procedures for doing fourth order factorizations. The fourth order convergence of the resulting algorithm allowed very large time steps to be used. In simulating the Brownian dynamics of 121 Yukawa particles in two dimensions, the converged result of a first order algorithm can be obtained by using time steps 50 times as large. To further demostrate the versatility of our method, we derive two new classes of fourth order algorithms for solving the simpler Kramers equation without requiring the derivative of the force. The convergence of many fourth order algorithms for solving this equation are compared. Typeset using REVTEX 1
منابع مشابه
0 Ju n 20 00 Higher Order Force Gradient Symplectic Algorithms
We show that a recently discovered fourth order symplectic algorithm, which requires one evaluation of force gradient in addition to three evaluations of the force, when iterated to higher order, yielded algorithms that are far superior to similarly iterated higher order algorithms based on the standard Forest-Ruth algorithm. We gauge the accuracy of each algorithm by comparing the step-size in...
متن کاملFourth-order algorithms for solving the multivariable Langevin equation and the Kramers equation.
We develop a fourth-order simulation algorithm for solving the stochastic Langevin equation. The method consists of identifying solvable operators in the Fokker-Planck equation, factorizing the evolution operator for small time steps to fourth order, and implementing the factorization process numerically. A key contribution of this paper is to show how certain double commutators in the factoriz...
متن کاملar X iv : m at h - ph / 0 00 70 23 v 1 1 5 Ju l 2 00 0 First order ODEs , Symmetries and Linear Transformations
An algorithm for solving first order ODEs, by systematically determining symmetries of the form [ξ = F (x), η = P (x) y+Q(x)] where ξ ∂/∂x+η ∂/∂y is the symmetry generator is presented. To these linear symmetries one can associate an ODE class which embraces all first order ODEs mappable into separable through linear transformations {t = f(x), u = p(x) y + q(x)}. This single ODE class includes ...
متن کاملSolving a Stochastic Cellular Manufacturing Model by Using Genetic Algorithms
This paper presents a mathematical model for designing cellular manufacturing systems (CMSs) solved by genetic algorithms. This model assumes a dynamic production, a stochastic demand, routing flexibility, and machine flexibility. CMS is an application of group technology (GT) for clustering parts and machines by means of their operational and / or apparent form similarity in different aspects ...
متن کاملar X iv : a st ro - p h / 00 06 29 0 v 1 2 1 Ju n 20 00 Gas Rich Dwarfs from the PSS - II III . H I Profiles and Dynamical Masses
We present Arecibo neutral hydrogen data on a sample of optically selected dwarf galaxies.
متن کامل